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Exact counting of Euler Tours for generalized series-parallel graphs 
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Abstract 

o 

f»o ■ We give a simple polynomial-time algorithm to exactly count the number of Euler Tours (ETs) of 

any Eulerian generalized series-parallel graph, and show how to adapt this algorithm to exactly sample 
a random ET of the given generalized series-parallel graph. Note that the class of generalized series- 
parallel graphs includes all outerplanar graphs. We can perform the counting in time 0(mA 3 ), where 
A is the maximum degree of the graph with m edges. We use 0(mA 2 log A) bits to store intermediate 
values during our computations. To date, these are the first known polynomial-time algorithms to count 
or sample ETs of any class of graphs; there are no other known polynomial-time algorithms to even 
approximately count or sample ETs of any other class of graphs. The problem of counting ETs is known 
to be (IP-complete for general graphs (Brightwell and Winkler, 2005 [3|) and also for planar graphs 

Q i (Creed, 2009 El). 
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1 Introduction 

Let G = (V, E) denote an undirected, connected multigraph where the degree, d{v), of each vertex v € V 
is even. Any standard introductory graph theory text has a result stating that every graph with even degree 
has an Euler tour (or Euler circuit), i.e., a circuit that traverses every edge of G exactly once. This result, 
which also holds for Eulerian multigraphs, implies a very simple linear-time algorithm for testing whether 
a given multigraph admits some Euler tour. In this paper we consider the counting and sampling of Euler 
tours of generalised series-parallel graphs, a special class of multigraphs. Throughout the paper we will use 
the term 'graph' to include multigraphs (ie, to allow the possibility of loops and parallel edges). 

In 2005, Brightwell and Winkler showed that the problem of counting Euler tours is (JP-completeEJ. 
This is in sharp contrast to the case of directed Eulerian graphs (i.e. connected digraphs for which the 
indegree equals the outdegree at each vertex), where the number of Euler tours can be counted exactly in 
polynomial-time using the Matrix-Tree Theorem [2] and the so-called "BEST" Theorem (after de Bruijn, 
van Aadenne-Ehrenfest, Smith, and Tutte, although apparently the first two deserve credit as the original 
discoverers HI). More recently, Creed |0]] showed that counting the number of Euler tours in undirected 
graphs remains (jP-complete if G is restricted to be a planar graph. 

These tJP-completeness results naturally lead one to question for which classes of graphs can (exact or 
approximate) counting of the number of Euler tours be done efficiently. We consider the case of generalized 
series-parallel (GSP) graphs, a subclass of planar graphs with distinguished source s and sink t vertices, that 
may be constructed in an inductive manner using a small number of operations. The key operations combine 
two generalized series-parallel graphs G\ , G2 to form a larger graph, and these operations are known as 
series composition G\o s G2, parallel composition G±o p G2, and the dangling composition G\OdG2- These 
operations are all defined in Section [2~T1 The problem of checking whether a given graph G is a generalized 
series-parallel graph can be done in polynomial time [6], and if G is a GSP graph, a hierarchical binary tree 
decomposition of G (see Section l2~TI ) can be found in polynomial time ||5][l2l[T3][l4l[T0l[6]. m a hierarchical 
tree decomposition T of a generalized series-parallel graph G, each internal node u is associated with an 
operation o u (one of o s , o p and Od) and each leaf node u is associated with an edge (of G). For every 
internal node u of the tree decomposition, with child nodes v and w, the subtree T u represents a connected 
subgraph G u of G, where o u is the top-level operator combining G v and G w . 

In this paper we give a polynomial-time "dynamic-programming"-like algorithm for exactly counting 
Euler tours for GSP graphs. Our main result is as follows: 

Theorem 1 Let G be an Eulerian generalized series-parallel graph having m edges and maximum degree 
A. We assume that we know T, a binary tree decomposition for G. 

Counting the number of Euler tours of G can be performed using C(mA 3 ) arithmetic operations, and 
using 0(mA 2 log A) bits for storing intermediate values in the computations. 

We also show how we can use our results to sample an Euler tour of a GSP graph exactly uniformly 
at random. To the best of our knowledge, these are the first polynomial-time algorithms for counting (and 
sampling uniformly at random from) the set of Euler tours for any class of undirected graphs. At the time 
of writing, we do not even know of any approximate-counting results for any significant class of undirected 
graphs. One of the few relevant results of which we are aware is an approximate asymptotic formula for the 
number of Euler tours of the complete graph on any odd number of vertices [ 8 ]. 

For general graphs, there exist natural Markov chains for sampling Euler tours such as the "Kotzig 
chain", which uses circuit reversals. However, although there has been some research on the Kotzig Markov 
chain [ 1 1], no correct proof of rapid mixing has yet been found for any class of Eulerian graphs. 

Our paper is structured as follows: in Section [2] we give some key definitions, including the definition 
of an (s, ^-decomposition of an Eulerian graph with distinguished vertices s,t. Most of the work in this 



paper deals with the relationship between these (s, ^-decompositions and Euler tours, and the building-up of 
counts of these (s, ^-decompositions. In Section[3]we will show that the number of Euler tours of a Eulerian 
graph G can be expressed as a simple weighted sum over the count of (s, ^-decompositions with k non-loop 
paths for a linear number of k values. Therefore the main component of both our counting and sampling 
algorithms will be to build a table containing the counts of (s, ^-decompositions, for all component graphs 
of G, and all relevant k. In Section|4]we show how the counts of (s, ^-decompositions of two GSP graphs 
G\, G2 (for varying values of k) can be used to build the counts for the parallel-composition G\o p G2, the 
series-composition G\o s G2 and the dangling-composition Gio^G^. Finally in Section [5] we sketch the 
simple polynomial-time algorithms that allow us to exactly count and sample Euler tours of GSP graphs. 

2 Definitions 

Throughout this paper we assume that G = (V, E) is a multigraph, which may contain parallel edges and 
loop edges. We will assume that the edges of the graph have some arbitrary but fixed ordering e±, . . . ,e m 
(where m is the number of edges of the graph). We will use Adjc{v) to denote the set of edges adjacent to 
the vertex v. We also use the notation No = {0, 1, 2, . . .}. 

2.1 Generalized series-parallel graphs 

First we define the class of graphs that we study in this paper. 

Definition 2.1 A generalized series-parallel graph (GSP graph) is any graph G = (V, E, s, t) with two 
distinguished nodes s,t G V, which can be built inductively from the following operations: 

B: The graph consisting of two vertices connected by a single edge is a GSP graph (where s and t are 
the endpoints of the single edge). 

o s : Given two GSP graphs G\ = (V\,Ei, s, t), G2 = (V2, E2, s', t'), the series composition ofG\, G2 is 
defined as G = Gio s G 2 =def (V, E, s, t'), where V = def (Vi U V 2 ) \ {>'} and E = def (E 1 UE 2 U 
{(t, v):Ve = (s', v) G Adjo 2 (s')}) \ Adj G2 (s'). 

o p : Given two GSP graphs G\ = (V\,Ei,s, t), G2 = (V2, E2, s', t'), the parallel composition ofG\, G2 
is defined as G = G\o p G2 = (V, E, s, t), where V =def (Vi U V2) \ {s', t'}, and E =def (E\ U^U 
{ef = (s,v) : Ve = (s',v) G Adj G2 (s')} U {e> = (t,v) : Ve = (t',v) G Adj G2 (t')}) \ {Adj G2 {s') U 
Adj G2 (t>)). 

Od-' Given two GSP graphs G\ = (Vi,Ei,s, t), G2 = (V2, E2, s', t'), the dangling composition ofG\, G2 
is defined as G = G 1 o d G 2 = (V,E,s,t), where V = de f (ViUV" 2 )\{«'}, and E = de f (EiUE 2 VJ{e' = 
(s,v):Ve = (s>,v)eAdj G2 ( S ')})\Adj G2 (s'). 

Intuitively, the series-composition of G\ and G2 is formed by identifying the sink t of Gi with the 
source s' of G2, with the new source and sink of G\o s G2 being, respectively, the source of G\ and the sink 
of G2. The parallel-composition of G\ and G2 is formed by identifying the source s of G\ with the source s' 
of G2, and the sink t of G\ with the sink t' of G2. The dangling-composition of G\ and G2 is formed by 
identifying the source s of Gi with the source s' of G2 (and keeping the source and sink of G\). Note that 
the order of G\, G2 is important for the o s and Od operations. Moreover, G\o s G2 and G20 S G\ need not be 
isomorphic to one another. 
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Figure 1: Examples of the operations for constructing GSP graphs. 



The class of series-parallel graphs consists of those that may be obtained using the three operations 
B,o s , and o p . Adding the dangling-composition operation brings us to the class of generalized series- 
parallel graphs. Out of interest we note that outerplanar graphs are known to be generalized series-parallel 
graphs [7]. Figure Q] shows some examples of the operations used to construct GSP graphs. 

We note here that GSP graphs have a compact representation using a binary tree. Each leaf of the binary 
tree is a (labelled) edge of the graph, and each internal vertex of the tree represents a series, parallel, or 
dangling operation that applies to its two children, and each internal vertex is labelled appropriately. In the 
case of a series (resp. dangling) operation, we can define the tree so that the left child corresponds to G\ and 
the right one G 2 in the definition of the series (resp. dangling) operation. See Figure [2] for an example. We 
haven't explicitly indicated the source and sink node for each of the subgraphs, but they are obvious for this 
small example. 

2.2 Euler tours and "legal" graphs 

Definition 2.2 A connected graph is said to be Eulerian if every vertex has even degree. A graph is said 
to be near-Eulerian (with discrepancies at v, v') if all but two vertices v, v' have even degree. We will say 
that a graph G = (V, E, s, t) with distinguished vertices s, t is legal if the graph is either Eulerian or is 
near-Eulerian with discrepancies at s and t. 

The particular class of legal graphs we will focus on in Sections 0] and [5] are GSP graphs. The four 
operations B, o s , o p and o^ for constructing GSP graphs only change the degree of the distinguished ver- 
tices s, t. Therefore, in building a GSP graph which is Eulerian, we need only consider graphs which are 
either Eulerian, or are near-Eulerian with discrepancies at s and t. Note that any graph constructed by the 
base-case operation B is by default a legal graph. A series-composed graph G = G\o s G 2 is a legal graph if 
and only the following conditions hold: 

si. G\ and G 2 are both legal graphs, and 

s2. either dc 1 (t),da 2 (s') are both even, or da r (t),da 2 (s') are both odd. 

A parallel-composed graph G = G\o p G 2 is a legal graph if and only if: 

pi. G\ and G 2 are both legal graphs. 
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Figure 2: A GSP graph construction in (a), and its corresponding binary tree decomposition in (b). 



A dangling-composed graph G = Gio^G^ is a legal graph if and only if: 

dl. G\ is a legal graph and G2 is an Eulerian legal graph. 
Note that for any legal graph, the parity of the source is the same as the parity of the sink. 

Observation 2 Let G = (V, E, s, t) be any legal GSP graph, and T any (rooted) binary tree decomposition 
of G according to o p , o s and 04- Then for every node u of T, the subgraph G u corresponding to the 
subtree T u is a legal GSP graph. 

Definition 2.3 For any Eulerian graph G, an Euler tour T is any path T in G which traverses every edge 
exactly once. We consider Euler tours to be identified under the operations of rotation and reversal. For any 
graph G, we let ET(G) denote the set of Euler Tours ofG. 

Definition 2.4 Let G = (V, E, s, t) be a legal graph, and let p be a (not necessarily simple) path in G. We 
say that p is a (s, t)-simple path if one endpoint of p is s and the other is t, but none of the intermediate 
points of p lie in {s,t}. An (s, s)-loop is a (not necessarily simple) circuit that starts and ends at s, such 
that neither s nor t are intermediate nodes of the circuit. A (t, t)-loop is defined similarly. 

Note that an (s, t)-simple path, an (s, s)-loop or a (t, t)-loop may visit vertices of V \ {s, t} more than once. 

Definition 2.5 Let G = (V, E, s, t) be a legal graph. An (s, t) -decomposition of G is any collection C = 
{pi, . . . ,pi} of edge-disjoint (s, t)-simple paths, (s, s)-loops, and (t, t)-loops in G such that U i=1 {e : e £ 
Pi} = E, and such that for every i G [£], the initial edge ofpi has an edge label with a lower index than that 
of the final edge ofp^ 



Informally, an (s, t) -decomposition is a partition of the edges of G into (s, t)-simple paths, and loops 
that contain either s or t, but not both. Any Euler tour of G gives rise to a unique (s, ^-decomposition in a 
natural way. Conversely, any fixed (s, ^-decomposition will give rise to a number of Euler tours in G (see 
Lemma[6]below). We have the following observation. 

Observation 3 Let G = (V, E, s, t) be a legal graph. If C is an (s, t)-decomposition of G, then £ = \C\ = 

d(s)+d(t) 
2 

3 (s, ^-decompositions and Euler Tours 

In this section we demonstrate the relationship between the set of (s, ^-decompositions of an Eulerian graph 
(with distinguished vertices s, t) and the set of Euler tours of that graph. 

Definition 3.1 Let G = (V, E,s,t) be an Eulerian graph with distinguished vertices s and t, and let T € 
ET{G). We say that the (s, i)-decomposition C ofG is consistent with the Euler tour T if for every p G C, 
either p or rev{p) is a contiguous segment ofT (where rev(p) is the reverse of the path p). 

Observation 4 Let G = (V, E, s, t) be an Eulerian graph with distinguished vertices s and t and suppose 
T £ ET(G). Then there is exactly one (s, t) -decomposition of G which is consistent with T. 

We make the following definition for all legal graphs. 
Definition 3.2 Let G = (V, E, s, t) be an legal graph. Let 

K (G) = {k G N | < k < mm{d(s), d{t)}, k mod 2 = d(s) mod 2}. 

Let C(G, k) denote the set of (s, ^-decompositions of G in which there are k (s, t)-simple paths. We define 
1 {G,k) = \C{G,k)\. 

Using simple counting and parity arguments, we have the following observation. 

Observation 5 For any legal graph G = (V, E, s, t), C(G, k) = %for any k k(G). 

Lemma 6 Let G = (V, E, s, t) be an Eulerian graph with distinguished vertices s, t, let k € k(G), and 
let C be an (s, t) decomposition of G such that C € C(G, k). Then the set of Euler tours of G which are 
consistent with C is a one-to-one correspondence with tuples of the form 

,{2 ) ...,fe} 1 r{Hl,.. ) ^} ) ^ + l,..M }i ^ iiG{1) _ 1} ^ j (1) 

such that 

• tv, t, a are permutations on k — 1, ^ s l~ and - L~~ elements respectively; 



x £ No ' is a sequence of ^ non-negative integers which sums to 



d(s)-k . 
2 ' 



y G Nq ' 2 is a sequence of -^ non-negative integers which sums to -M; — . 



Proof: Suppose that C consists of k (s, t)-simple paths, k s (s, s)-loops, and kt (£, £)-loops. Then since the 
paths of C contain exactly one copy of each edge of G, it must be the case that k s = (d(s) — k)/2 and 
kt = (d(t) — k)/2. Moreover, k must be an even number. 

Let Pi,---,Pk be the (s,£)-simple paths in C; Pk+i, ■ ■ ■ ,P(k+d(s))/2 be the (s, s) -loops in C, and 
P(fc+d(s))/2+i, • • • ,P(d(s)+d(t))/2 be the (t,t)-loops in C. Assume wlog that p 1 is the path in C whose s- 
adjacent edge has the lowest edge label, of any of the paths p\, . . . ,pk- We will assume wlog that p\ is the 
initial path in every Euler tour, and that it is oriented from s to t, thereby enforcing the rule that an Euler 
tour is not changed by a rotation or a reversal. 

We will characterize all Euler tours consistent with C by initially considering the order in which P2, ■ ■ ■ ,Pk 
appear in the tour after p\ is traversed. This order is given by a permutation n on k — 1 elements. Note that all 
of the paths p k +i, ■ ■ ■ ,P(d(s)+d(t))/2 are (s, s)-loops or (t, t) -loops. Therefore if the ordering p^ 2 ),- • • ,Pw(k) 
is to be extendible to an Euler tour, then for every 2 < i < k, we must direct the path p n t i \ from s to t if and 
only if i is odd (and from t to s if i is even). Therefore it determines both the order and the direction of all 
the paths p2 , . . . , Pk in the Euler tours for that ir. 

To complete the Euler tour, consider all ways of inserting the loops Pk+i, ■ ■ ■ ,P(d(s)+d(t))/2 i nto the 
partial tour pi,p w (2), • • • ,Pn(k)- in what follows, assume that 7r(l) = 1. Also identify ir(k + 1) with 7r(l) 
and 7r(0) with ir(k). The (s,s)-loops can only be inserted into the intervals after some path p n ni) an< i 
before p 7r (2j+i), for some 1 < i < |. This implies that there are k/2 intervals where the (d(s) — k)/2 
(s, s)-loops can be inserted. We may insert as many as or (d(s) — k)/2 of the (s, s)-loops into any of 
these k/2 positions, in any order. Moreover, each of the (s, s)-loops may be oriented in either of the two 
possible directions. These choices may be encoded in terms of 

• A permutation r on the paths Pk+i, ■ ■ ■ , P d{s)+k , specifying the order of insertion. 

2 

• A sequence x G N of non-negative integers which sums to ^ s '~ , where X{ specifies the number 
of (s, s)-loops to be inserted between p^^i) an(1 PiT(2i+i)- 

d(s)-k 

• A vector b 6 {— 1, 1} 2 , specifying a direction for each of the (s, s)-loops. 

Similarly, the (t, t)-loops can only be inserted into the intervals after some path p n (2i-i) an d before p^^i) . 
for some 1 < i < |. There are (d(t) — k)/2 (t, t)-loops, and fc/2-positions where they may be inserted. 
These choices may be encoded in terms of 

• A permutation a on the paths p d( s )+k , , , . . . , p d( s )+d(t) , specifying the order of insertion. 

2 + 1 2 

k/2 d(t)~k 

• A sequence y € N of non-negative integers which sums to -Mj — , where m specifies the number 
of (t, t)-loops to be inserted between p^^i-i) an d P^{2i)- 

d(t)-k 

• A vector b €{—1,1} 2 , specifying a direction for each of the (t,t) -loops. 

We write b = b'b" to obtain a sequence of length (d(s) + d{t))/2 — k over {—1, 1}. 

Finally note that for any pair of Euler tours T, T' that are both consistent with the decomposition C, 
either T ~ T' under rotations and reversal (and therefore T and T' are the same tour), or else the tu- 
ple (ir(T),T(T),a(T),x(T),y(T),b(T)) induced by T will differ from the tuple (Tr(T'),T(T'),a(T'), 
x(T'),y(T'), b(T')) for T". This proves the one-to-one correspondence between the Euler tours with de- 
composition C and the tuples (ir,T,a,x,y,b). ■ 



Corollary 7 Let G = (V, E, s, t) be an Eulerian graph with distinguished vertices s, t. Then, for k(G) as 
given in Definition [ 



\ET(G)\ = £ (t)^2^)^))/^(^-l)!(^-l)! 7 (G^). (2) 

Proof: The number of k for which C(G, k) ^ is finite. Any C € C(G, k) must contain (d(s) — k)/2 (s, s)- 
loops and (d(t) — k)/2 (£, £)-loops in order to satisfy the conditions of an (s, ^-decomposition. Therefore 
we must have k < d(s) and k < d(t). Also, we require k mod 2 = d(s) mod 2 = 0. 

By Observation 0J every T € ET(G) is consistent with a unique (s, t) -decomposition of G. We can 
therefore express \ET(G)\ as a sum over the set of (s, ^-decompositions, of the number of Euler tours 
consistent with that decomposition. 

In Lemma [6] we proved that for any k € k(G), the set of Euler tours consistent with any (s,t)- 
decomposition C G C(G, k) is in one-to-one correspondence with the set of tuples described in (Q~|). The 
number of these tuples depends on k, on (d(s) — k)/2 and on (d(t) — fe)/2 but not on the specific decompo- 
sition itself. Therefore, we can express \ET(G)\ as a weighted sum of the j(G, k) values, where the weight 
against j(G, k) is the number of tuples for k with respect to G. There are (k — 1)1 possible it permutations, 
((d(s) — k)/2)\ possible r permutations, and ((d(t) — k)/2)\ possible a permutations. Counting all the 
possible x vectors is equivalent to counting the number of sequences of k/2 non-negative numbers which 
sum to (d(s) — k)/2. This can be expressed as the number of ways of partitioning (d(s) — k)/2 into k/2 
parts, which is ( (d(s) ~S/^I ( i /2)_1 ) = ^ffcMwf 1 ) ■ Similarly the number of y vectors satisfying the speci- 
fied constraints is (' k/L^i )• ^ ne numDer of 6-vectors is 2( d ( s ) +(i ( t ))/ 2_fe . Therefore the number of Euler 
tours which are consistent with C, for any C € C(G, k), is exactly 



2 J V 2 

_i) 
(fc-l)! x 2 



(^1- 1)1(^1-1)1 
\ 1)I(|-1)! 



*.■ 



^'(T-V) 


/'-M 1 — A d(s) + d(t) 

(I-J 2 ! 


2 2 ft . 





Then using the identity (&) | = rTT^fpjVp » we § et tne desired expression for \ET(G)\. ■ 

Corollary |7] gives an explicit formula to evaluate the number of Euler tours, given the value of j(G, k) 
for every feasible k. In the next section we will show how to recursively compute the j(G, k) values of a 
legal generalized series-parallel graph. 

4 Recursive computation of (s, ^-decompositions 

In Section[3l we have an exact formula for the number of Euler tours of a Eulerian graph, given the j(G, k) 
values for all feasible k. We now show how the j(G, k) values for a legal GSP graph G can be recursively 
computed using the GSP structure of G. In Section [5] we apply these recurrences to exactly count and 
exactly-sample Euler tours. 

4.1 Parallel combination 

Lemma 8 Let G\ = (V\,Ei,s,t) and G2 = (V2,E2,s',t') be two legal GSP graphs and consider the 
parallel composition G\o p G2, with source s ~ s' and sink t ~ t'. Let k € k(Gio p G2)- The (s,t)- 
decompositions of G\o p G2 with C € C{G\o p G2-,k) are in one-to-one correspondence with pairs of the 



form (Cx,C 2 ), where C\ G C(Gx,kx), C 2 G C(G2,k — ki) for some k\ such that kx G k[G\) and 
k-fae k(G 2 ). 

Proof: Let C G C(Gxo p G2,k). By construction of G±o p G2, we know that any (s,t)-simple path, (s,s)- 
loop or (t,t)-loop of G\o p G2 either lies entirely in Gx or entirely in G2. Define C\ to be the set of paths 
and loops of C which lie entirely in G\ . Define C2 to be the set of paths and loops (with s' substituted for s) 
which lie entirely in G2. The pair (Cx, C 2 ) is unique by this definition. Then C\ is a (s, t) -decomposition 
for G\ for some k\ G k{G\), k\ < k, and C2 is an (s 1 , t') -decomposition with k — k\ (s' , t')-simple paths. 

On the other hand, consider any pair {C\, C2) where C\ G C(G\, kx) and C2 G C(C2, k — kx), such that 
k\ G k{Gx) and k — kx G ^(G^). Then C\ U C2 G C{G\o p G2, k\ + fo), where C2 is a copy of C2 with 
s', t' replaced by s, t respectively. 

These two parts prove the one-to-one correspondence. ■ 

Corollary 9 Let G\ = (Vi, Ex, s, t) and G2 = (V2, E2, s' , t') be two legal GSP graphs and consider the 
parallel composition G\o p G2, with source s ~ s' and sink t ~ t'. Let k G k{Gxo p G2)- Then the number of 
(s, ^-decompositions ofGxo p G2 which contain k (s, t) -simple paths is 

j(Gxo p G 2 ,k)= Y, l(Gx,kx)*7(G 2 ,k-kx). 

ki€K,(Gi),k-ki€K(G2) 

4.2 Series combination 

The recursive relationship between the set of (s, t')-decompositions of G\o s G2 and the (s, t) -decompositions 
of Gx and (s', t') -decompositions of G2 is more involved than for the parallel case. We will see in LemmafTOl 
that the decompositions which are used to form (s, t') -decompositions of C(Gxo s G2, k) are the elements of 
C(G t , kx) for kx > k and of C{G 2 , k 2 ) for k 2 > k. 

Definition 4.1 Suppose we are given two series parallel graphs Gx = (Vx ,Ex,s,t) and G2 = (V 2 ,E 2 ,s',t') 
such that ddis) mod 2 = dG 2 { s ') mod 2. Suppose kx G k(Gx) and k2 G k(G 2 ). We define D(kx,k 2 ) 
(with respect to Gx, G 2 ) to be 

D{kx,k 2 ) = d Gl (t) + d G2 {s') - kx - k 2 . 

Lemma 10 Let Gx = (Vx, Ex,s,t) and G 2 = (V 2 , E 2 , s', t') be two legal GSP graphs such that d Gl (s) mod 
2 = dc 2 (s') mod 2. Consider the series composition Gxo s G2 of the two graphs, with source s and sink t'. 
Then for any k G n(Gxo s G2), there is a one-to-many correspondence between the elements ofC(Gxo s G2, k) 
and between tuples of the form 

U U (d. C ^ *iU • • • ^1 ^{1, . . . , k 2 }, x G N ( fcl+fe)/2 , a{l, ..., ^%^}, {1, -l}^.fe)/2 

fe 1 6K(G 1 )fc 2 6«(G 2 ) 
ki>k k2>k 

such that 

• Cx is an (s, t)-decomposition from C(Gx, kx) and C2 is an (s', t')-decomposition from C(G 2 , k 2 ). 

• 7Ti, 7T2, <r are permutations on kx elements, /^ elements and - j 2 ' elements respectively. 

• x is a non-negative integer vector of length (kx + k 2 )/2 which sums to D(kx, k 2 )/2. 
The factor in the one-to-many correspondence is $ s (fe, kx, k 2 ) =d e f kl^f^\^ z 2 z ^\2 ( ' kl+k2 ^ 2 ^ k . 



Proof: First consider an arbitrary (s, t') -decomposition C of G\o s G2 with C G C(G\o s G 2: k). Recall that 
t ~ s' is the point where G\ was joined to G2 ■ For every path p G C, let pi , . . . , p^( p ) be the set of segments 
of p such that for each pi, the endpoints, but none of the intermediate points of pi lie in {s, t', t ~ s'}, 
and such that each segment is directed such that the initial edge has the lower edge label. Then for every 
path p G C and every 1 < i < £(p), we know that pi either lies entirely in G\ or entirely in G 2 . We define 

d = Uj, 6C {Pi:l<J<^),PiGGi} 
C 2 = U pe cfe : 1 < i < i(p),Pi G G 2 }. 

Ci is then an (s, ^-decomposition of G\ and C2 is an (s', £') -decomposition of G 2 . Suppose that C was an 
element of C(G\o s G 2 , k) for k G k(Gio s G2). Then each of the k (s, i')-simple paths of C will contribute 
one (s, t)-simple path to C\ and one (V , t')-simple path to C 2 . This is because every path between s and t' 
in GiOgG^ must pass through t ~ s'. Hence if Ci G C(G±, fa) and C2 G C(G 2 , k 2 ), we are guaranteed that 
k < k\ and k < k 2 . By our rule for directing the pi segments, the pair (C±, C 2 ) is uniquely defined for any 
given C. 

Now suppose we are given C\ G C(G\,fa), C 2 G C(G 2 ,k 2 ) for fa G n{Gi),k 2 G k(G2). We will 
characterize the (s,t') -decompositions of G\o s G 2 which can be constructed from C\ and C 2 . First of all 
note that by the argument of the previous paragraph we can only use (Ci, C 2 ) to construct decompositions 
which contain at most min{fci, k 2 } paths with both s and t' as endpoints. Suppose we are considering such 
a k G k,{G\o s G 2 ), k < min{A;i, k 2 }. 

In Ci we have k\ (s,t)-simple paths, {dc^s) — k\)/2 (s,s)-loops, and {dc^t) — &i)/2 (t, i)-loops. 
In C 2 we have k 2 (s', t')-simple paths, (dG 2 (s') — k 2 )/2 (s', s')-loops, and (dG 2 (t') — k 2 )/2 (t',t') -loops. 
For C, we need to construct k (s, t')-simple paths, (dc 1 (s) — k)/2 (s, s)-loops, and {dc 2 (t') — k)/2 (f, t')- 
loops. We make the following observations: 

(s, t')-simple paths: By construction of G\o s G 2 , every path between s and t' in G±o s G 2 must pass through 
t ~ s'. Therefore, every (s, t')-simple path in C must be built using exactly one of the k\ (s, t)-simple 
paths of C\, and exactly one of the k 2 (s', t')-simple paths of C 2 . 

(s, s)-loops: Each of the (s, s)-loops of C\ must become an (s, s)-loop in C. The extra number of (s, s)- 
loops that we will need to add to C is (dci i s ) ~ k — (dd (s) — k\))/2 = (k\ — k)/2. 

(t', t')-loops: Each of the (£', t')-loops of C 2 must become a (t f , t')-loop in C. The extra number of (£', t')- 
loops that we will need to add to C is (d G2 (t') - k - (d G . 2 (t') - k 2 ))/2 = (k 2 - k)/2. 

Observe that after k paths have been chosen from the (s, t)-simple paths of C\ and from the (V, t')-simple 
paths of C 2 , the paths and loops of C\ U C 2 which have not been allocated any role in C are as follows: 

• k\ — k remaining (s, t)-simple paths from C\, 

• k 2 — k remaining (s 1 , t')-simple paths from C 2 , 

• all of the (d Gl (t) - fa)/ 2 (t, t)-loops from C\, 

• all of the (d G2 (s') - k 2 )/2 (a', s')-loops from C 2 . 

Neither the (t, t)-loops from C\ nor the (V, s')-loops from C 2 contain either of the distinguished vertices s 
and t' of G\o s G 2 ; therefore they play no significant role in constructing the extra (s, s)-loops nor the extra 
(£', t')-loops needed for C. The (fa — k) remaining (s, t)-simple paths from C\ must therefore be used 
to construct the extra (fa — k)/2 (s, s)-loops for C. Hence, we only need to construct a pairing of these 
remaining (s, t)-simple paths that were not used for the (s, t')-simple paths of C. Similarly, in any C which 



is built from the paths in C\ U C 2 , we must use the remaining (k 2 — k) (s', t')-simple paths to construct the 
extra [k 2 — k)/2 (t', t')-loops needed for C. Again, any pairing of these (k 2 — fc)-paths is sufficient. 

We now observe that the choice-and-pairing of the k (s,t)-simple paths from C\ and the k (s',t')- 
simple paths from C 2 to form the (s, t')-simple paths of C can actually be combined with the pairing of the 
remaining {k\ — k) (s, t)-simple paths, and also the pairing of the remaining {k 2 — k) (V, t')-simple paths, 
as follows: we permute all of the (s, t)-simple paths of C\ (a permutation wi of length ki) and also permute 
all of the (V, t')-simple paths of C2 (a permutation of length k 2 ). Then we construct the following pairings: 

(s,t') -pairing: We pair the 7Ti(i)-th (s,t)-simple path of C\ with the 7T2(i)-fh (V, t')-simple path of C2, 
for every 1 < i < k. 

(s, s)-pairings: We already have ((^(s) — k)/2 (s, s)-loops from C\. The extra (k\ — k)/2 (s, s)-loops 
are constructed by pairing ni(i) with m(i + 1) for k < i < k\, i mod 2 = (k + 1) mod 2. 

(£', £')-pairings: We already have (d G2 (t') - k)/2 (f , t') -loops from C 2 . The extra (k 2 - fe)/2 (£', f)-loops 
are constructed by pairing ^(i) with ^(i + 1) for k < i < k 2 , i mod 2 = (k + 1) mod 2. 

Observe that in this model, the number of ways of coming up with the same set of pairings is k\ ( -^— ) ! ( -^— ) ! 
2( fe i+ fe 2)/2-fc observe that for all these pairings, our requirement to order the paths of C in terms of lower- 
edge first implies that the choice of the pairing determines the relative order of the pairs, for every pairing 
we have constructed. 

Finally, the (t, t)-loops of C\ and (s', s')-loops of C 2 must be included in our decomposition C. There 
are D{k\, k 2 )/2 of these in total. These paths may be inserted at any point where t appears in the partial 
paths we have constructed so far. There are exactly (ki + k 2 )/2 occurrences of t ~ s' in the (s, t'), (s, s) and 
(t, £)-pairings we have constructed at this point. We may insert as many as or all D(k\, k 2 )/2 of the (t, t) 
and (s 1 , s') loops into any individual position - the particular partition chosen is encoded as the x vector. 
Each of the D(ki, k 2 )/2 items is different, so once the partitioning has been determined, there are 2 2 \ 
ways of ordering the (i, t)- and (V , s')-loops for insertion - encoded by the a permutation. Finally, no (£, £)- 
or (s', s')-loop q will ever be inserted so that it is adjacent to an endpoint of any path of C. Therefore for 
every such q, the 2 directions of inserting q (q and rev(q)) result in a different C-decomposition. ■ 

Corollary 11 Let G\ = (Vi, E\, s, t) and G 2 = (V 2 ,E 2 , s', t') be two legal GSP graphs such that d Gl (s) mod 
2 = da 2 (s') mod 2. Consider the series composition G\o s G 2 of the two graphs, with source s and sink t'. 
Then for any k £ k(Gio s G 2 ), the number of (s, ^-decompositions ofG\o s G 2 with k {s,t} -simple paths is 

7 (G l0s G 2 ,fc) 



kl lk 2 l ( ^-^ -l)! 2 ( dGl (Wc 2 (s'W 



D E 7(G r i,fc 1 )*7(G2,*2)* iblMlMr (h+lti x, >- 2kl+k2 . k 



fcl6«(Gi)fc 2 6«(G 2 ) • 2 • 2 • I 2 

fcl>fc k 2 >k V 

Proof: We use LemmaQIJl For any particular k\,k 2 , there are fci! m permutations, k 2 \ tt 2 permutations, and 
{D{k\ , k 2 )/2)l a permutations. The number of x vectors is the same as the number of ways of partitioning 
D(k\, k 2 ) items into {k\ + k 2 )/2 parts, which is 

(k 1 + k 2 )/2 + D( kl ,k 2 )/2-l\ _ ( dGl(f)+ 2 ^ (s ° -l)! 



(k 1 + k 2 )/2-l )~ ( ki+k 2 1 ^ g(fc 1 ,fc 2 ) ,- 

The number of tuples for a particular C\,C 2 with C\ € C(G\,ki) and C 2 G C(G 2 , k 2 ), is therefore 



h\k 2 \ 



( 2 -1)! D{h,k 2 ) ^ D(klM)/2 



•fei+M _ i\] D ( k iM) \ 
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Dividing by k\^\^\^^-\2^ kl+k2 ^ 2 ~ k , and cancelling some terms, we find that the number of (s,t') 
decompositions of G\o s G2 with k (s, t')-simple paths that can be constructed from C\, G 2 is 

hlh[ ^ Gl (t)+d G2 (s') _ 1) , 2(dGi(f ) +dga (, 0)/2 



k\( kl ~ k )\( fc 2~ fc V / fci+fc2 _ ]\j 2 fc i+ fc 2-fe 



4.3 Dangling combination 

The dangling combination changes the degree of the source and leaves the degree of the sink unaffected in 
a GSP graph. In terms of the (s, ^-decomposition, the dangling operation adds (s, s)-loops at the source s. 

Lemma 12 Let G\ = (V±,Ei, s, t) be a legal GSP graph and G 2 = (V2, E2, s' , t') be an Eulerian GSP 
graph and consider the dangling composition Gio<iG 2 , with source s ~ s' and sink t. For every k G 
n{G\OdG2) \ k(Gi), C(GiOdG2,k) = 0. For every k G k(C?i), the elements of C(GiOdG2,k) are in 
one-to-many correspondence with tuples of the form 

|J (c lj C W l,..., te }, xeN ^V{l,...,*^ 

fc 2 GK(G 2 ) 

such that 

• C\ is an (s, t) -decomposition from C(G±, k) and C 2 is an (s', t') -decomposition from C(G 2 , /c 2 ) where 
k 2 G k(G 2 ). 

• 7r 2 , a are permutations of length /c 2 and {dG 2 {t') — &2)/2 respectively. 

• x is a non-negative integer vector of length fc 2 /2 which sums to (dc 2 (t') ~ & 2 )/2. 

The factor in the one-to-many correspondence is < J ) d(A; 2 ) =def -jf ! 2 2 ' 2 . 

Proof: First observe that for any k, and any C G C(Gio^G 2 , A;), C induces an (s, ^-decomposition C" on 
Gi with exactly fc (s, t)-simple paths. Therefore we must have C(GiOrfG 2 , fc) = for every k $ k(G\). 

Now, for any /c G «(Gi) consider an arbitrary (s, ^-decomposition Gof GiOrfG 2 withG G C(GiOrfG 2 , fe). 
Recall that s ~ s' is the point where Gi is joined to G 2 . For every path p G G, let pi, . . . ,p« p ) be the 
set of segments of p such that for each pi, the endpoints, but none of the intermediate points of pi lie in 
{s(~ s'),t, t'}, and such that each segment is directed such that the initial edge has the lower edge label. 
Then for every path p G C and every 1 < i < £(p), we know that pi either lies entirely in Gi or entirely 
in G 2 . We define 

Ci = U pe c{Pi ■ 1 < i < £(p),Pi G Gi} 

C 2 = U p6C { K :l<i<f(j)),K£G 2 }. 

Gi is then an (s, t) -decomposition of Gi and G 2 is an (s', t') -decomposition of G 2 . Suppose that G is an 
element of C(Gio,iG 2 , k) for A; G k(Gio^G 2 ). Then a path p G G is an (s, t)-simple path of G if and only 
if it is an (s, t)-simple path of C\. This is because every (s, £)-simple path of G has to lie entirely in G\ and 
hence is an (s, t) simple path in C\. This, in turn, implies that C\ G C(G\, k) and G 2 G C(G2, fc 2 ) where fc 2 
is independent of A; and fc 2 G k(G 2 ). By our rule for directing the pi segments, the pair [C\, G 2 ) is uniquely 
defined for any given G. 

Now suppose we are given C\ G C{G\, k\) and G 2 G C(G2, fc 2 ) where fei G rc(Gi), /c 2 G rc(G 2 ). We 
will characterize the (s, t) -decompositions of Gio^G 2 which can be constructed from Gi and G 2 . First of 
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all note that by the argument of the previous paragraph, an (s, ^-decomposition of G±OdG 2 constructed 
using (Ci, C 2 ) will have exactly k\ (s, t)-simple paths. In other words, if C G C{G\OdG 2 , k) is constructed 
using (Ci, C 2 ), then C\ G C(Gi, &i) where k\ G k(Gi) and fci = k. We also note that a path p G C is a 
(£, t)-loop in C if and only if p is an (i, t)-loop in C\. This is because the vertex t lies in V{G\) \ V(G 2 ). 

In C\ we have fci (s,£)-simple paths, (dc 1 (t) — k\)/2 (t,t)-loops, and (cfe 1 (s) — fei)/2 (s,s)-loops. 
In C2 we have fo (s', t')-simple paths, (dG 2 (t r ) — k 2 )/2 (£', £')~1°°P S > an d (^g 2 ( s ') ~~ &2)/2 (s', s')-loops. 
For C, we need to construct fc(= k\) (s, £)-simple paths, (dc 1 (t)—k)/2 (t, £)-loops, and (dc 1 (s)+dc 2 i s ') — 
k)/2 (s, s)-loops. We make the following observation: 

(s, s)-loops: Every (s, s)-loop of C\ and every (s', s')-loop of C 2 must become an (s, s)-loop in C. The 
extra number of (s, s)-loops that we will need to add to C is {do 1 (s) + dc 2 ( s ') — k)/2 — (do 1 (s) — 
k)/2-(d G2 (s')-k 2 )/2 = k 2 /2. 

Observe that the (s, t)-simple paths and (t, t)-loops of C\ have become the (s, t)-simple paths and (t, t)- 
loops of C and the (s, s)-loops of C\ and C 2 have become the (s, s)-loops of C. The paths and loops of 
C\ U C 2 which have not been allocated any role in C are as follows: 

• all of the &2 {s', t')-simple paths from C2, 

• all of the (d G , 2 (t f ) - k 2 )/2 (f, t')Toops from C 2 , 

All of the k 2 (s',t') -simple paths and the (dG 2 (t') — k 2 )/2 (t 1 , t') -loops from C 2 would be combined to 
form A;2/2 (s, s)-loops in C ( extra (V, s')-loops in C 2 ). As the (t' ', t')-loops do not contain the vertex s', we 
ignore them for the time being. To obtain the ^2/2 (V, s')-loops, we pair up the k 2 (s' , t')-simple paths in 
the following fashion. We permute all the (s' , t')-simple paths (a permutation ir 2 of length k 2 ) and construct 
the following pairing: 

(s', s')-pairing: The extra (V, s') loops are constructed by pairing the 7T2(i)-th (s', t')-simple path of C 2 
with the ir 2 (i + l)-th (s 1 , t')-simple path of C 2 , where 1 < i < k 2 and i = 1( mod 2). 

Observe that in this model, the number of ways of coming up with the same set of pairings is -^\ 2 k2 ' 2 . 
Observe that for all these pairings, our requirement to order the paths of C in terms of lower-edge first 
implies that the choice of the pairing determines the relative order of the pairs, for every pairing we have 
constructed. 

Finally, the (£', t')-loops of C 2 must be included in our decomposition C. There are (da 2 (t') — k 2 )/2 
of these in total. These loops may be inserted at any point where t' appears in the partial (s', s')-loops we 
have constructed so far. There are exactly ^2/2 occurrences of t' in the extra ^2/2 (V, s')-pairings we have 
constructed at this point. We may insert as many as or all {dG 2 (t') — k 2 )/2 of the (t' ', t') loops into any 
individual position - the particular partition chosen is encoded as the x vector. Each of the (dc 2 (£') — k 2 )/2 

items is different, so once the partitioning has been determined, there are 2 ~ ' wavs °f ordering the 

(t', t')-loops for insertion - encoded by the a permutation. Finally, no (£', i')-loop q will ever be inserted so 
that it is adjacent to an endpoint of any path of C. Therefore for every such q, the 2 directions of inserting q 
(q and rev(q)) result in a different C-decomposition. ■ 

Corollary 13 Let G± = (V\,Ei, s, t) be a legal GSP graph and G 2 = (V 2 ,E 2 , s', t') be an Eulerian GSP 
graph. Consider the dangling composition G\OdG 2 of the two graphs, with source s and sink t. Then for 
any k G K,{G\OdG 2 ) = k{G\), the number of (s,t) -decompositions of G\OdG 2 with k {s,t} -simple paths 
isj(Gio d G 2 ,k) 

-7(0.,*). E «'««s#ii*(¥ 1 - 1 ) ! * 2d ' ,gW/I " fe 
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Proof: We use Lemma [T2l For any given fo there are k 2 \ ^2 permutations and ((d(5 2 (t') — &2)/2)! cr 
permutations. The number of 3; vectors is the same as the number of ways of partitioning (do 2 it') — k 2 )/2 
items into k 2 /2 parts, which is 

k 2 /2 + (d G2 (t>) - k 2 )/2 - 1\ _ (^ - 1)! 



(fe/2) - 1 / " (*a _ ^1 / ^(t')-fc2 \ , • 



2 •*•/• v, 2 
The number of tuples for aparticular C\, C 2 with Ci € C(G\, fc) and C2 € C(G 2 , k 2 ), is therefore 



'/fca _ 1 nj / d G Jt')-fc 2 \, 



Dividing by -^\ 2 k2 ' 2 , and cancelling some terms, we find that the number of (s, ^-decompositions of 
G\OdG 2 with k (s, £)-simple paths that can be constructed from C±, C 2 is 

k 2 \ fd G2 (t' 



/»G 2 VM _ 1 \ U2 d G2 (t')/2-k 2 



£!(£-!)! V 2 



5 Algorithms 

Our counting and sampling algorithms will use a compact binary tree representation T(G) for the GSP 
graph G, as in Figure |2l This representation can be computed in polynomial-time (see [6]), and is of size 
linear in the number of edges of G. 

5.1 Counting Euler tours 

Consider a binary tree representation T of a given Eulerian GSP graph G = (V, E, s, t). We assume that 
every vertex u € V(T) represents a subgraph G u of the graph G, obtained by applications of o s ,o p , or o^ 
to the graphs within the subtree at rooted at u. As noted in Observation |2j such component graph G u in a 
hierarchical decomposition of G will satisfy the property of being a legal GSP graph. For each vertex u in T 
let s u denote the source of G u and t u the sink of G u . 

For each such u G V(T), we compute the values of j(G u ,k) for k G k(G u ). To do this, for each 
u G V(T), we define variables d u (s u ), d u (t u ) and the integer array 7«(fc), for 1 < k < m (where m = \E\) 
with each vertex u of the tree. The value d u (s u ) will be used to store the degree of the source s u in G u and 
the value d u (t u ) will be used to store the degree of the sink t u in G u . For every k G k(G u ), j u {k) will be 
used to store the value of j(G u , k). All values for the internal nodes are initialized to 0. We compute the 
values for d u (s u ), d u (t u ) and j u (k) for all nodes u G V(T) in a bottom-up fashion. 

Each leaf of T(G) corresponds to a graph consisting of a single edge of the graph G, ie, to a B operation. 
Hence we set d u (s u ) = 1 and d u (t u ) = 1 for every leaf u. For a graph G u consisting of a single edge, we 
have k(G u ) = {1}. Therefore we set 7 U (1) = 1 and J u (k) = for k ^ 1. 

Alternatively, if u is not a leaf of T, then G u was obtained from the two smaller legal GSP graphs 
corresponding to the child nodes of u via the o s or the o p or the Od operation. Let I and r denote the left 
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and right children of u. We set the values of d u (s u ) and d u (t u ) as follows: 

Series operation : s u = sg ,d u (s u ) = d^(s^) , t u = t r ,d u (t u ) = d r (t r ) 
Parallel operation : s u = sg = s r , d u (s u ) = d^(sp) + d r (s r ) , 
t u = t(. = t r , d u (t u ) = dg(tg) + d r (t r ) 
Dangling operation : s u = sg = s r ,d u (s u ) = d^{sp) + d r (s r ) , 
t u = te,d u (t u ) = di(ti) 

In the case of o p , we use the equation of Corollary [9] to compute the value of 7 u (/c) from the Jg(-) and 7 r (-) 
values, for every k G k(G u )- In the case of o s , we use the equation of Corollary QT] to compute the value 
of 7«(fc) from the Jg(-) and 7 r (-) values. In the case of o^, we use the equation of Corollary [T3l to compute 
the value of 7 u (/c) from the 7^(-) and 7 r (■) values. 

Finally, once we have the values of j(G, k) for k G n(G) for the given GSP graph G, we can use 
equation © to compute \ET(G)\. 

5.2 Time and space complexity of counting 

Given an Eulerian GSP, G, what is the time complexity of computing \ET(G)\1 We first assume that we 
already have a binary tree decomposition that describes how to construct G (as described in Section 12.11) 
using the standard operations for building GSPs. Note, first, that if |£?(G)| = m, then the binary tree has m 
leaves and m — 1 internal vertices, as each operation in the construction of G combines two legal subgraphs 
to make one new connected (legal) subgraph. 

Corollary |7] tells us how to find \ET(G)\, given the values of j(G,k) for all k G n(G). Let A = 
max degree of G. From Corollary [71 we see that the number of terms in the sum for computing \ET(G) \ is 
at most A/2 (as k must agree with the parity of the degree of the source s). So we need to know how much 
time is required to compute the j(G, k) values for any fixed k G k{G). We consider each operation in turn, 
and examine the time needed to compute the values 7(Gio x G , 2, k) for an operation o x G {o p , o s , o^}, given 
all the corresponding values for the two graphs G\ and G2. 

Parallel combination: Forafixedfc G K(Gio p G2),fromCorollary|9]we see there are min{|«;(G ! i)|, |«(G ! 2)|} 
terms in the sum to find r y{G\o p G2, k), given the values for each of G\ and G2. Therefore, computing 
j(Gio p G2,k) takes 0(A) operations. (Note: Here we are assuming that we before beginning our com- 
putations to find \ET(G)\, we first perform a "pre-processing step" by computing, and storing, values k\ 
and 2 k for < k < A. We can do this "pre-processing" in time 0(A), using 0(A 2 log A) and 0(A 2 ), 
respectively, bits of space for the k\ and 2 k values.) 

Series combination: Corollary ITTItells us there are at most \k(G-\ )| ■ |«;(CM| G 0(A 2 ) terms in the (double) 
summation for finding 7(Gio s C?2), k for a fixed k G k{G\o s G2)- 

Dangling combination: Finally, from Corollary [13] we see there are |k(G2)| G 0(A) terms in the summa- 
tion for finding 7(GiOrfG2, k) for a fixed k G ^(GiOdG^). 

Overall, we see that computing any value j(H, k), for any legal subgraph H of G (corresponding to an 
internal vertex of the binary tree decomposition), and any fixed k G k(H), takes time 0(A 2 ). Since there 
are at most 0(A) values of k G k(H), and m — 1 internal vertices in the binary tree decomposition, we 
can compute all j(H, k) values with 0(mA 3 ) operations. Finally, since there are 0(A) terms in the sum 
to find \ET(G)\, we see we can find the number of Euler tours of a GSP in time 0(toA 3 ). 
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Taking a crude upper bound of O ( ( 2 ) ( 2 ) • • • ( 2 ) j G (A A ) for the number of Euler tours of an 

Eulerian graph with max degree A, we see that we need at most 0( A log A) bits to store one j(H, k) value. 
Using, again, that there are 0(A) values of j(H, k) at each of the m — 1 internal vertices of the binary tree 
decomposition of G, we find we need 0(mA 2 log A) bits to store all values necessary to compute |£T(G)|. 
The bounds above on the time and space complexity of computing \ET(G)\ are those given in Theo- 
rem [TJ 

5.3 Sampling Euler tours 

We now show how to sample uniformly from ET(G). We will sample a Euler tour by first sampling an 
(s, t) -decomposition, for some k € j(G, k), and then applying Lemma[6]to obtain a uniform random tour. 
Note that after counting the number of Euler tours of G, we now know value of^ u (k) = j(G u , k) for every 
component graph in the tree decomposition. 

(s, ^-decompositions are grouped in terms of the number of (s, i)-simple paths, k, in the decomposition. 
Therefore we first must choose k € k(G), with the appropriate probability. Using Corollary 13 it is clear 
that the probability that an element of ET(G) is consistent with an (s, t) -decomposition with k (s, t)-simple 
paths, for any given k € k(G), is exactly 

il)p- k l(G,k) 

E Iei .( G )(i)^-' 7 (G,o' 

Our sampling algorithm computes these probabilities for all k € k(G), and then chooses k € k(G) to 
be an exact sample from this known distribution. 

Having fixed k for G, we must choose k\ and k2 values for the graphs G\ and G 2 that combine (either 
via o s or o p or o^) to form G. Suppose G = Gio s G 2 . We require (k\,k 2 ) sucn that k\,k 2 > k. The 
probability with which we pick any (&i, k 2 ) such that ki € K,(Gi) and ki > k for i = 1, 2 is exactly 



'<j(ti)+d(g 2 ) 



7(Gi,fci)*7(G 2 ,A;2)* 



1 ! 



k^.k^. \ I "'j' 2 (d(t 1 )+d( S2 ))/ 

..i fc l- fc i fc 2- fc ; 7 k 1 + k 2 ^T ' 2 k 1 +k 2 -k 

7(GTfc) 



2 



For the case of o p , the probability of choosing any (fci, k — k\) for k\ € k,{G\), k — k\ € ^(^2) is 

7(Gi,fci)*7(G2,fc-fci) 

7 (G,fc) 

For the case of Od, the probability of choosing any (k, k 2 ) for k G k(G\), k 2 G ^(^2) is 

7(Gi,fc)*7(G 2 ,feO* fc9 t' *(^-l>^2 (t ' )/2 -^ 
-TfK-^-- 1 ) 1 v / 



(3) 



(4) 



7(G,fe) 



(5) 



Having chosen ki, we recurse in the respective sub-trees and use ([3]>, (01) or © to choose k values with the 
appropriate probability for every internal node in the binary decomposition tree. Then for every u S T(G), 
we have generated a value k u from the exact uniform distribution for the (s u ,t u ) -decomposition of G u . 
Next we must generate the (s u , t u ) -decompositions themselves. 

Suppose at node u we have G u = G\o p G 2 and we have C{ the respective (sj, tj) -decompositions for 
graph Gi, i € {1, 2}. Note that s\ ~ s 2 and t\ ~ t 2 . We need to construct the (s u , t u ) -decomposition for 
G u . From Lemma[8l it is clear that {C\, C 2 ) is an (s u , t u ) -decomposition for G u . 

Suppose at node u we have G u = G\o s G 2 . Let d be the (sj, tj)-decompositions generated for graph 
Gi, i G {1, 2}, where C\ is a exact uniform sample from C(Gi, ki) for the fixed fej values, i = 1,2. Note 
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that t\ ~ S2,s M ~ s± and t u ~ £2- Then by Lemma [TOl we can combine C\ and C2 in a suitable fashion 
to obtain an (s u , £ M ) -decomposition for G u . We have to ensure that every (s u ,t u )- decomposition with k 
(s u ,t u ) simple paths which have both s u and t u as end-points that could arise from the combination of the 
given decompositions of G\ and G 2 is equally likely. By Lemma [TOl we know that we can construct a 
random (s u , t u ) -decomposition of G u from C\, C 2 by generating a random tuple of the following form: 

m{l, ..., h}, n 2 {l, ..., k 2 },x G N f 1+fc2)/2 , a{l, . . . , ^%M } , {1> _ l} D(k lM )/2 

such that the vector x sums to D(k\, k2)/2, and then following the steps described in the proof of LemmafTOl 
This random tuple can easily be generated in polynomial-time. Hence we have the o s case. 

The (s u , t u ) -decomposition in the dangling case can also be constructed by first recursing on the (smaller) 
generalized series-parallel graph in a similar fashion as above. (Recall that we know that this smaller graph 
is itself Eulerian.) We use this decomposition of the subgraph to form a set of (s u , s u ) loops that we add to 
the decomposition for the parent node (in the tree) using the method described in Lemma ITS] 

The last and final step is to combine the (s, t) -decomposition at the root node of the binary decompo- 
sition tree into an Euler tour where s and t are the terminals of the generalized series-parallel graph. This 
can again be done simply, and in polynomial-time, by generating a random tuple of the form described in 
Lemma [6] and following the steps in the proof of that Lemma. 

5.4 Other Counting Problems in GSP graphs 

We wish to point out that the following list of combinatorial structures could be counted exactly using the 
above technique:(i) independent sets, (ii) matchings, (iii) A;-colourings (A;-constant), (iv) dominating sets. At 
the same time, we do not lay claim to being the first to do so. In fact, a recent result by Steve Noble [9] 
shows that (i),(ii) and (iii) could be counted exactly for a larger class of graphs, namely, the class of bounded 
treewidth graphs. 
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